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1. Introduction 

Supersymmetry (SUSY), is a well motivated Beyond the Standard Model (BSM) scenario. 
Electroweak precision data indicates that SUSY should be realized at the TeV scale or 
below. If this is the case, it will be accessible to direct experimental measurements at 
the Large Hadron Collider (LHC) through the production of SUSY particles. Monte Carlo 
simulations have shown the possibility of the discovery of TeV-scale SUSY with 1 fb -1 of in- 
tegrated luminosity [1,2]. First measurements involving supersymmetry-sensitive variables 
have been performed by both ATLAS [3-6] and CMS [7] collaborations. 

At the LHC, colored particles like the SUSY partners of quarks and gluons, i.e. squarks 
and gluinos, will be copiously produced. Theoretically, these processes are extensively stud- 
ied within the Minimal Super symmetric extension of the Standard Model (MSSM). The 
leading-order (LO) contributions, of 0(a s ), have been known for a long time [8-12]. The 
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next-to leading order (NLO) QCD corrections have been computed [13, 14] and imple- 
mented in the public code PROSPIMO [15]. They affect the LO predictions substantially 
and they reduce the scale dependence considerably. More recent is the estimation of 
the logarithmically enhanced next-to-next-to-leading order (NNLO) QCD contributions 
to squark hadroproduction, the resummation of the QCD Sudakov logarithms at the next- 
to-leading-logarithmic (NLL) accuracy, and the resummation of the leading Coulomb cor- 
rections [16-22]. Their contribution amounts up to 10% for squark and gluino masses of 
the order of 1 TeV. They further stabilize the prediction against scale variation. 

The electroweak (EW) contributions have been computed for several processes pro- 
ducing colored SUSY particles. They exhibit an extremely rich and complicate pattern, 
in particular when the EW contributions appear already at tree level. The tree-level EW 
contributions are of 0(a 2 + a s a). At the parton level, they arise from gg-annihilation, 
(/(/-scattering or photon-induced processes, depending on the final state considered. They 
are known [23-28] and they can increase the LO cross section by up to 20%. The im- 
pact of these contributions in the context of non-minimal flavor violation and explicit CP 
violation has been investigated as well [29,30]. The NLO EW corrections contribute at 
0{a 2 s a) and have been computed for stop-anti-stop [23,31], squark-(anti-)squark [28,32], 
gluino-squark [25], and gluino-gluino [33] production. Their size is comparable with that 
of the tree-level EW and NNLO QCD contributions, and their impact strongly depend on 
the SUSY scenario considered. 

The production of third generation squarks is special. The non-negligible mixing in 
the stop and sbottom sector could lead to relatively low masses for the lightest bottom 
and top squarks, favoring their direct production at the LHC. Moreover 6-tagging makes 
bottom- and top-squark production experimentally distinguishable from the production of 
the squarks of the first two generations [34-36]. A dedicated analysis looking for third 
generation squark production at the LHC is already available [6]. This kind of searches are 
particularly important in SUSY scenarios such as the ATLAS benchmark scenario SU6, 
where inclusive searches with jets, missing transverse energy, and leptons are problem- 
atic [2]. 

In this paper we focus on the hadronic production of bottom-squark pairs 

PP^b a b},b a b p ,b* a b}, a,Pe{l,2}. (1.1) 

In particular we present the first complete computation of the NLO EW corrections to 
diagonal sbottom-anti-sbottom pair production, 

PP^b a b* a , a €{1,2}. (1.2) 

The contribution of the remaining processes is small (cf. Section 5), hence we will not 
include them in our discussion on the EW corrections. The process ( |1.2j ) exhibits specific 
features like the mixing between left- and right-handed b-squarks, the renormalization of 
the sbottom sector [37,38], the non-negligible Higgs-boson contributions, and the enhanced 
Yukawa couplings for large values of tan f3 with the related need of resummation [39] . These 
features make the computations of the electroweak contributions to the processes ( |1.2j ) 
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substantially different from those for squark-anti-squark [32] and stop-anti-stop [23] pro- 
duction, and justify a specific investigation, which is reported in this paper. 

The outline of the paper is as follows. In Section |2| we summarize the tree- level contribu- 
tions to the processes ( |1.2j ). Section [| describes the various partonic processes contributing 
at 0(a 2 a) and the strategy of the calculation. The numerical impact of the NLO EW 
contributions at the LHC with \/S = 14 TeV and \/S = 7 TeV is presented in Section |j. In 
Section [| we discuss the numerical impact of the subleading bottom-squark pair production 
processes. The Feynman diagrams and the technical details of the renormalization of the 
sbottom sector are collected in the Appendix. 



2. Tree-level cross section 



In this section we describe the tree- level contributions to the process (|1.2j), which are 
of order 0(af), O (a s a), and 0(a 2 ). We will conventionally denote the cross section 
(amplitude) of a partonic process X at a given order 0(a°a fe ) as da a / (MY). The 
parton luminosities are defined as 



X 



(2.11 



where ft(x,^p) is the parton distribution function (PDF) of the parton i inside the 
hadron A. 

2.1 Tree-level QCD contributions 

The leading-order cross section, of the order 0(a 2 ), is given by 

-1 at rl 



<K>° Q ? £ (S) = [ dr ^da 2 >° --(§) + V [ dr ^da 2 ' - - (s) 



T ~' — " a J T 

f j dL b i 2,0 
J dr bb^b a b% y 



(2.2) 



where tq = 4m~ /S is the production threshold. S and s = tS are the squared center-of- 

b a 

mass (cm.) energies of the hadronic and partonic processes, respectively. The sum runs 
over q = u, d, c, s. The three classes of partonic processes contributing are 

g(pi)g(p2) K(p3)K(p4), (2.3a) 
q(pi)q(p2)^b a (p 3 )b* a ( P4 ), (2.3b) 

b(Pi)HP2)^b a ( P3 )b* a (p 4 ). (2.3c) 
The corresponding partonic cross section can be obtained from the Feynman diagrams in 



Figure E^l. In terms of the Mandelstam variables, 

s = {pi+P2? 1 i=( Pl -p 3 ) 2 , u=( yPl - Pi ) 2 , (2.4) 
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the differential partonic cross section for a given subprocess ££' — > b a b* a can be written as 

with the squared lowest order matrix element averaged (summed) over initial (final) state 
spin and color. 



2.2 Tree-level EW contributions 



The tree-level electroweak (EW) contributions, which are of the order O(a s o) and 0(a 2 ), 
read as follows, 



PP->b a b; 



(S) 



dr^Ida 1 ' 1 - - (s 

dr 91^b a b* a ■ 



dL 



+ / dr 



dL 



(2.6) 



bb 



TO 



dr 



d^ 1 - r,(s)+daj 2 - ~(s) 

bb-+b a b± v y bb^b a b% y ' 



The contributions of O(o 2 ) arise from the processes ( |2.3b ) and ( 2.3c ). The partonic cross 
sections, 



da ' 2 



qq-^b a l 



{s) = Y\m 



0, 1 

qq^iab"^ 



2 dt 



16lTS 2 



do"?r 2 r r„(s) 

bb^b a b* y ' 



El" 



0,1 

bb^tb a b* a 



2 dt 



167TS 2 



(2.7) 



are obtained from the diagrams in Figure B.l. In the case of process (2.3c), the diagrams 
with t-channel gluino and neutralino exchange further allow for a non-vanishing QCD-EW 
interference term of 0(a s a), 



da 1 : 1 - . (s) = 2YRe\M 1 : - - (m : 1 - - )*} 

bb^b a b* a K > v bb^b a b* a v bb^b a b* a j / 



dt 



167TS Z 



(2- 



We approximate the CKM matrix V with the unity matrix, thus we consistently neglect 
the diagram depicted in Figure plKd). The latter dia gram, together with the diagrams in 
Figure B.l| (b), gives rise to 0(a s a + a 2 ) contributions. These contributions are at least 



quadratic in \V c f,\ ~ 10 • \ V u b\ ~ 4 • 10 and numerically negligible (see Section |4.3| ) 



The first 0(a s a) contribution in Eq. (2.6) arises from the photon-gluon induced pro- 



cess 



g{Pi)l(P2) 



(2- 



The corresponding partonic cross section can be obtained from the diagrams in Figure B.2 
and reads as follows 



da 



1. 1 



,(*) = E|" 



1/2,1/2 



2 dt 



167TS 2 



(2.10) 
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3. Next-to-leading order EW contributions 

In this section we list the NLO EW corrections to the process (|1.2| ). These contributions 
are of the order 0(a 2 a) and arise from virtual corrections and bremsstrahlung processes. 
Using an obvious notation, the corresponding contributions to the total cross section read 
as follows 



d(J NLO EW 



PP^b a b, 



AS) 



dr 



dL 



</</ 



dr 



j.2,1 , 
gg->bj>* 



da*' 1 - r Js) + dd'- - T ,_ [8] 



r 2 ' 1 - - ( 

gg^b a b* a ~i 



+E / drdL "~ 



n Jt 



dr 
dL, 



dr gq->b a b*q 



.„(*) + 



dL, 



5 9 _ d ~2,i 



dr m-^bab^q 



We do not consider the contributions arising from the bremsstrahlung processes 
70i) q(P2) -> b a (p 3 ) b* a {p 4 ) q(p 5 ), j{pi) q(p 2 ) -> b a (p 3 ) b* a (p 4 ) q(p 5 )- 



(3.1) 



(3.2) 



As already pointed out in Ref. [32], they are suppressed because of the 0[a) suppression 
of the photon PDF inside the proton. Moreover, these processes are further suppressed 
by an additional factor a s with respect to the process ( |2.9[ ) and thus negligible. The 
0{a 2 s a) contributions of the partonic processes with a bottom quark in the initial state are 
neglected as well. The reason is twofold. First of all these contributions are suppressed by 
the bottom PDF with respect to the contributions in Eq. (|3.1|). In addition they have an 
additional factor a s with respect to the 0(a s a + a 2 ) contributions of the process ( |2.3d ), 



which turn out to be small (cf. Section 4.3). 

The amplitudes are generated and algebraically simplified with support of FeynArts [40, 
41] and FormCalc [41,42], while the numerical evaluation of the one- loop integrals has been 
performed using LoopTools [42]. Infrared (IR) singularities are regularized giving a small 
mass A 7 and X g to the photon and to the gluon, respectively. The mass of the light quarks 
is kept in order to regularize the collinear singularities. 

3.1 Virtual corrections 

The 0(a 2 a) contributions to the partonic process ( |2.3aD are given by 



d^ 1 -- (s) = 2YRe\M 1 >° - - (Mh\. )*} 
gg^b a b* a y > L^i y gg^b a b* a \ gg-+bj>* a ) J 



di 



167TS 2 



(3.3) 



where A4 1,0 is the tree- level amplitude while M 1,1 is the one- loop amplitude obtained from 
the diagrams depicted in Figure |B.3| . The virtual corrections to the process (|2~3b| ) read as 
follows 

^- b 5. (*) = *J>lM\° r - v (M 1 , 1 - - V +M '} - - (m\° - - V) 

qq^b a b^ s y qq^b a b* a \ qq^tb a b* a J qq^b a b* a \ qq^b a b* a J J 167TS 2 

(3.4) 
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A4 0,1 and .A/f 1,0 are the tree-level EW and the tree-level QCD amplitudes, respectively. 
A4 1 ' 1 is the one-loop amplitude obtained from the EW insertions to the leading-order 
diagrams and from the QCD corrections to the tree- level EW diagrams (Figure |B,4| ). .A/f 2 ' 
is the amplitude corresponding to the QCD box diagrams depicted in Figure [B.5| . 

In order to cancel the UV divergences we need the 0(a) renormalization of the wave- 
function of the light quarks and of the sbottom sector. The field renormalization constants 
are fixed in the on-shell scheme, in analogy to that described in Ref. [32]. The renormal- 
ization of the sbottom sector has to be performed together with that of the stop sector. 
In order to avoid numerical instabilities and artificially big contributions from the coun- 
terterms, care has to be taken in choosing the renormalization scheme [37,38]. We use the 
"DR bottom-quark mass" scheme introduced in Ref. [37]. Since in particular regions of the 
MSSM parameter space this scheme can give rise to numerical instabilities, we have explic- 
itly checked its reliability in the SUSY scenarios considered in this paper (cf. Section |4|). 
The explicit expression of the renormalization constants in the "DR bottom-quark mass" 
scheme are collected in Appendix |A[ 

3.2 Real corrections 

The O{a 2 o) contributions to the partonic processes with a photon in the final state, 

g{Pi) g{P2) -> K(p 3 ) b* a (p 4 ) 7(p 5 ), (3.5a) 
q{Pi) q(P2) -> ba(ps) KiPi) 7(P5)j (3.5b) 



are obtained from the tree-level diagrams in Figure |B.q . The phase space integration 
is divergent in the soft photon region, i.e. in the region p® — > 0. In the case of the 
process flOED further singularities arise in the collinear region, i.e. p\^, • P5 — > 0. IR and 
collinear singularities are treated using the phase-space slicing method. The description of 
the method and the relevant formulae are collected in Ref. [28]. 
The gluon bremsstrahlung process, 

q(Pl) qip%) ->• b a {p 3 ) b* a (p 4 ) g(p 5 ), (3.6) 

contributes at 0{a 2 a) via the interference of QCD-based and EW-based Feynman diagrams 
depicted in Figure [B.7| . IR singularities of gluonic origin are treated in close analogy to the 
photonic case. Color correlations are taken into account using the formulae collected in 
Appendix B of Ref. [28]. Due to the color structure, the interference term of a QCD-based 
and an EW-based diagram vanishes if both gluons are emitted from an initial-state or a 
final-state particle. 

Real quark radiation contributes at 0(a 2 a) as well, 

g(pi) q{P2) ->■ b a (p 3 ) b* a (p 4 ) q(p 5 ), (3.7a) 
g(pi) q{P2) ->■ b a (p 3 ) b* a (p 4 ) q(p 5 ). (3.7b) 

This IR- and collinear-finite set is given by the interference of QCD and EW tree-level 



diagrams (cf. Figure B^). Only the interference from initial-state and final-state radiation 
contributes. 
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The IR singularities arising in the gg channel cancel in the sum of virtual corrections, 
process ( |3.3|) , and real photon radiation ( |3.5a ). In the qq channel the sum of the virtual 
corrections (^OJ) and of the contributions of real photon radiation ( |3.5b| ) and real gluon 
radiation fl3,6| ) is IR finite. This sum is affected by universal collinear singularities of 
photonic origin that can be absorbed in the PDFs. This can be achieved by means of the 
following substitution [43], 



fq{x,H F ) -> Ux,H F ) 1 



ae: 



7T 



cte q 



l-Ss 



dz 

z 



(3.8) 



e„ is the electric charge of quark q expressed in units of the positron charge, while 



1 - In 5 S - In 2 5. + In S s + - In 



a4 

m 2 



1 



2vr 2 



+ 7 9+ — + 31n^-21n z 5 S , 



K c (z) = P qq {z) In 



l4 L_ 

m 2 (1 — z)' 



P qq {z) In 



1 - z 



3 1 



21-2 



+ 2z + 3 



(3.9) 



with the splitting function P qq {z) = (1 + .z 2 )/(l — z). The factorization is performed in 
the DIS scheme. The replacement of the PDFs in Eq. ( |2.1| ) gives further contributions 
of 0(a 2 a) to the total cross section. As already mentioned, they cancel the collinear 



singularities affecting the O(a 2 s o) contributions (3.1). 



3.3 Resummation in the b/b sector 

The Higgs sector in the MSSM corresponds to a type-II two-Higgs doublet model, i.e. 
the down- type quarks couple to H\ and the up- type quarks to Hi- After spontaneous 
symmetry breaking, the up- (down-)type quarks get their mass by their coupling to the 
vacuum expectation value V2 (v\) of H2 (Hi). At tree- level, the bottom-quark mass m& is 
related to the H\bb Yukawa coupling via 

mb = \bV\. (3.10) 

Radiative corrections induces and effective H^bb coupling that can significantly alter the 
tree- level relation ( 3.1 0| ) [39,44-47]. This higher-order contributions do not decouple at 



low energies and are enhanced by a factor tan (3 = v 2/v 1 . As shown in Ref . [39] , the leading 
tan/3 enhanced terms can be resummed by using an appropriate effective bottom-quark 
Yukawa coupling. We follow Ref. [37] and use an effective Yukawa coupling defined as 
follows, 

1 mffa) + m b Am b _ m™^ 

A 6 = -1 , A = > (-a- 11 ; 

v\ 1 + Amj, v\ 
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where Amj is given by 



2a 



A, 2 



167T 



r ^/M2 tan /3 



cos 2 %/(m fl ,M 2l |i) + sin 2 9 h I{m h , M 2 , //) 



1 



1 



+ ^ cos 2 ^I(m Si , M 2 , /x) + - sin 2 ^ I(m^ , M 2 , /i) 



I(a,&, c) 



1 

( a 2 _^)(62 _ c 2)( a 2 _ c 2) 



(3.12) 



n 2 h 2 r 2 

2l2i a , 1.2 2 1 v , 22i c 

a b log—T + bc log ^ + c a log ^ 

o z c z ar 



Large logarithms from the running of the Yukawa coupling A at the renormalization scale 
jiR are resummed using the DR bottom-quark mass, 



m^ s is the on-shell bottom-quark mass defined according to 



(3.13) 



m ps = m M»( mz )6f 



MS/ 



ishift 



(,shift 



1 + 5. 



MS\2 



(3.14) 



S fin ' is the finite part of the scalar self-energies defined according to the Lorentz decom- 



position ( |A.7 ) of the Appendix [A|. The term proportional to Amj in the numerator of 



Eq. ( 3.11 ) has to be inserted to avoid double counting of the one- loop contributions of the 
resummed terms. 

Further tan f3 enhancement effects arise from three-point functions involving Higgs- 
bottom vertices. The tan /3-enhanced terms can be taken into account by modifying the 
Tibb coupling 5^5. 1 The combined effect of the resummation in the relation between A and 
m and of the resummation in the Higgs-bottom vertices is accounted for by performing 
the following substitutions, 



9h%b -> 9h°bb\ x ,^\ b I 1 



Ann, 



9hH 



9HHb\ Xb ^- Xb I 1 + Am 6 



tan /3 tan a 
tana 



tan /3 J ' 



9G°bb ->■ 9c°bb ■ 



t .a„/3V' (3 ' 15a) 
(3.15b) 



The coupling involving the Goldstone boson G° is not modified since the contribution 
from the vertex corrections exactly compensates the contribution of the bottom- Yukawa 
coupling resummation. 

4. Numerical results 

In this section we perform a detailed numerical analysis for diagonal sbottom-pair produc- 
tion at NLO EW. We stick to the notation introduced in Ref. [28]. The leading order cross 



"H stands for any of the neutral Higgs and Goldstone bosons, i.e. H — h°, H°, A° 7 G°. 
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section, the tree-level EW and the NLO EW contributions to the cross section are labeled 
by 

cJ Born = CJ 2,0 ; A(J tree EW = (fJ l, 1 + ^0, 2 )? A(J NLO EW = ^2, 1 ^ 

respectively. Acr EW = Acr tree EW + A<r NLO EW will be referred to as the EW contribution. 
The total sum of the LO cross section with the EW contributions is denoted by <7 NLO = 
<r Bom + A<7 EW . Relative EW contributions are defined by 

jtree EW _ Acr trcc EW ^Born ^NLO EW _ A(J NLO EW /^-Born ^EW _ ^^EW / a Boln 

(4.2) 

In distributions 5 denotes the relative EW contribution defined as 5 = (Onlo - 0Bom)/0Born> 
where O is a generic observable and Onlo is the sum of the Born and the EW contributions. 

4.1 Input parameters 

The Standard Model input parameters are chosen in correspondence with [48,49], 

M z = 91.1876 GeV, M w = 80.399 GeV, 

a' 1 = 137.036, a s (M z ) = 0.119, (4.3) 

m t = 170.9 GeV, mf(m z ) = 2.94 GeV. 

The strong coupling constant a s has been defined in the MS scheme using the two-loop 
renormalization group equation with five active flavors. 

For the numerical analysis we consider the mSUGRA scenarios SPSla' and SPS4. The 
first one is a "typical" SUSY scenario proposed by the SPA convention for comparison 
with other calculations [50]. The scenario SPS4 is characterized by a large value of tan/3. 
Within this scenario we study the dependence of the total cross section on the squark 
masses and on tan/3. The third scenario considered is the GMSB scenario SPS8. 

The particle spectrum is determined following the procedure described in Ref. [28]. 
Starting from GUT-scale parameters, cf. Table [j], we use the program Softsusy [51] to 
evolve the soft-breaking parameters down to the SUSY scale Mstjsy • In accordance to the 
SPA convention, a common SUSY scale Msusy = 1 TeV has been chosen. To get the right 
mixing in the sbottom sector, we first compute the low energy masses and use these to 
calculate the effective bottom-quark mass m^ R ' eff , Eq. ( 3.11 ). This mass is then used in 



the bottom-squark mass matrix to calculate the sbottom-mass eigenstates. The lighter of 
the two bottom-squarks is taken as the dependent squark. Its mass is therefore fixed by 
SU(2) invariance. Table |2] collects the shift Am&, cf. Eq. ( 3.12| ) and the effective bottom- 



quark mass together with the on-shell mass of the bottom squarks, of the gluino and of the 
lightest neutralino and chargino. 

Unless otherwise stated, the results presented in this section are computed setting the 
hadronic center of mass energy to \fS = 14 TeV and using the MRST2004QED parton dis- 
tribution functions [52]. The factorization and renormalization scales are set to a common 
value, fjL = = = m l a -> i- e - to the mass of the produced bottom squark. 
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m 


mi/2 


A 


tan/3 


sign(/i) 


SPSla' 


70 GeV 


250 GeV 


-300 GeV 


10 


+ 


SPS4 


400 GeV 


300 GeV 





49.4 


+ 




A 




N 

x T mess 


tan (3 


sign(/x) 


SPS8 


100 TeV 


200 TeV 


1 


15 


+ 



Table 1: High-energy input parameters for the different SUSY scenarios considered. The mass 
parameters m , m 1 / 2 and A are given at the GUT scale, tan/3 is evaluated at Msusy = ITeV. 





Am/, 


DR,cff 

m b 


bi 


t>2 


9 




xt 


SPSla' 


0.037 


2.38 


500 


533 


609 


101 


180 


SPS4 


0.23 


2.05 


428 


633 


736 


123 


217 


SPS8 


0.03 


2.42 


1070 


1085 


141 


253 





Table 2: The shift Ami, and the resulting effective bottom-quark mass as well as the on-shell 
masses of the bottom squarks, of the gluino, and of the lightest ncutralino and chargino within 
the different SUSY scenarios considered. Amj and m^ R ' cS are evaluated at the scale used for b±bl 
production. All masses are given in GeV. 

4.2 Total hadronic cross section 

Table || shows the hadronic cross section for diagonal bottom-squark production within the 
three considered scenarios for y/S = U TeV. 

As expected, the total cross section is dominated by the LO QCD contribution of 
O{ol). The contribution of the photon- induced process is independent of the mixing 
angle. In each scenario, its yield relative to the leading order cross section is similar 
for the two processes considered. Although formally suppressed by a factor a s , the NLO 
EW corrections are typically bigger than the tree-level EW contributions. In the SPSla' 
(SPS4) scenario the tree level and NLO EW contributions are more important in case of 
bib* (62^2) production. This can be explained by the chirality dependence of the SU(2) 
coupling and by the fact that in the SPSla' (SPS4) scenario bi (62) is mostly left-handed. 
In the SPS8 scenario the bottom squarks are twice as heavy as in the aforementioned 
scenarios, thus the Born cross section is about two orders of magnitude smaller. Further, 
the mixing between left- and right-handed squarks is more important and the sbottom 
masses are nearly degenerate. These features partially soften the differences among the 
tree-level EW contributions to bib* production and the ones to 62^2 production. 2 Huge 
cancellations between the qq and the gg channel amplify the dependence of the NLO EW 
contribution on the production process considered. As a result the NLO EW contributions 
to bib\ production and the ones to 62^2 production have opposite sign, the latter being 
three times bigger than the former. Summing up the various contributions, the relative 
yield in the scenario considered is below 2%. 

Table |4| collects the hadronic cross section for \fS = 7 TeV. The leading order total 

2 In the no- mixing limit the tree-level EW contributions to 6l production is one order of magnitude 
larger than the one contributing to ba production. 
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Table 3: Hadronic cross section for diagonal b a b* a production at the 14 TeV LHC within three 
different scenarios. Shown are the LO cross section, the tree-level EW as well as NLO EW contri- 
butions and the relative corrections as defined in the text. The numbers in brackets refer to the 
integration uncertainty in the last digit and are omitted if the uncertainty is at least one order of 
magnitude smaller than the quoted precision. All cross sections are given in femtobarn (fb). 



cross sections are reduced proportionally to the mass of the produced squark. They amount 
to 1 - 10% of their value at y/S = 14 TeV. In all the scenarios considered, the contribution 
of the photon-induced channel is enhanced with respect to the 14 TeV case. In the SPSla' 
(SPS4) scenario the impact of the NLO EW contributions to the production of the mostly 
left-handed bottom-squark, b\b\ (62^2)1 i s reduced. In contrast, the NLO EW contributions 
become more important in case of the production of the mostly right-handed sbottom. 
In the SPS8 scenario the EW contributions of the various channels are enhanced. In 
particular the NLO EW contributions are positive for both production processes. This 
is a consequence of the enhancement of the NLO EW contributions to the gg channel at 
y/S = 7 TeV (cf. Section [D]). 

4.3 Parameter scan 

The impact of tan (3 and of the sbottom masses have been studied performing a parameter 
scan on these parameters. In this scan, the soft breaking parameters Ml, M~ h , and Mj 
appearing the squared mass matrix, Eq. (|A.2| ), are set to a common value w. S q Uar k. All 
other parameters are set to their SPS4 values. The scans presented in this section are 
obtained for three different values of f7i S q Uar k = {300,600,900} GeV. tan/3 is varied from 
10 to 50. 
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Table 4: Same as Table || but considering diagonal b a b* a production at the 7 TeV LHC. 



In all scenarios considered, we have verified the smallness of the bottom-initiated tree- 
level contributions, justifying our procedure of neglecting the O(a^a) contributions to this 
channel. In Figure [l] we show the relevance of the various production channels at tree- 
level. For b\b\ and 62^2 production the gluon fusion channel contributes 70 — 90% of the 
total cross section. The remaining 10 — 30% of the total cross section is given by the qq 
channel, its relative yield increasing with tan /3. The contribution of the bb channel is at 
the permille level, in accordance with the analysis in Ref. [21]. 3 Owing to the small yield 
of the bb channel at tree-level, we will safely neglect the NLO EW contributions to this 
channel. 

The reliability of the renormalization scheme in the scenario considered has been veri- 
fied by investigating the behavior of the dependent parameters since these can potentially 
get large finite shifts by their renormalization constants. In the case of m~ and we have 
checked that the finite part of their renormalization constant is smaller than the parameter 
itself. The mixing angle renormalization constant enters only the counterterm of the 6J62 
and b\ &i quadratic terms via the combination 

5Y b = (ml-ml)se b . (4.4) 

3 The big contributions from the bb channel quoted in Ref. [27] are a consequence of two enhancement 
factors. First of all, the Higgs exchanged in the s-channel is resonant. Moreover, the Hbb Yukawa couplings 
are enhanced by the choice of a negative value for the parameter n, such that Amj, « —0.76. In our analysis 
the Higgs bosons are not resonant and we do not consider negative values of /i, which are disfavored by the 
measured value of the anomalous magnetic moment of the muon (g — 2) M [53]. 
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Figure 1: Relative yield of the various tree- level production channels as a function of the common 
squark mass breaking parameter m squar k- The left (right) plot shows b\b\ (^2^2) production. 

0.4 



0.3 
0.2 
0.1 


-0.1 
-0.2 




600 



rl 

r2 
r3 



900 



rl 
r2 
r3 



10 



20 30 

tan/3 



40 50 10 20 30 40 50 10 20 30 40 

tan /3 tan (3 



50 



Figure 2: The ratios r\,T2, and r% as defined in Eq. (4.5) as a function of tan/3 for different values 
of the common squark mass breaking parameter m = TO S quark in GeV. 



We have checked that the finite part of <5Yf, is smaller than the mass of the bottom squarks. 
Figure shows the ratios 
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m~ 



(4.5) 



61 



as a function of tan j3 for the various values of w-squark- 

The value of r\, r2, and r% is below 

0.2 in all the scenarios. 

We have explicitly checked the impact of the 0(a s a + a 2 ) contributions arising from 
the diagram depicted in Figure |B.l| (d). In the region of the parameter space considered 
in this paper these contributions are negligible. Indeed its numerical value is six orders of 
magnitudes smaller than the Born cross section, justifying the approximation of identifying 
the CKM matrix with the unit matrix. 

The results of the scan are collected in Figure [3|. The qq contribution is the sum of 
the tree-level EW and of the NLO EW contributions from the qq annihilation channel. 
The peaks in the corrections correspond to neutralino, chargino, or sfermion thresholds. 
These unphysical singularities affect the self-energy of the produced sbottom and can be 
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Figure 3: Relative yield of the EW contribution for b\b\ production (left panels) and for 62^2 
production (right panels) as a function of tan/3. Shown are the various production channels as well 
as the combinded effect. The soft breaking parameter TO squar k is defined according to Section 4.3 
and is expressed in GeV. tan/3 is varied between 10 and 50. 
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regularized by taking into account the finite widths of the unstable particle [54]. The 
curves in Figure || exhibit a step-function like behavior in the region where the mass of the 
produced sbottom is above threshold. Both, for b\b\ and for 62^2 production the behavior 
of the various contributions strongly depend on the size of tan (3. In the following we 
will distinguish between the low and the high tan /3 region, which are separated by the 
threshold. 

The case of b\b\ production is shown in the left panels of Figure H The EW contribu- 
tions of the qq and the gg channel have substantially the same tan/3 dependence close to 
threshold. In the low tan /3 region the contributions from the gj channel cancel against the 
one coming from the qq annihilation. The overall effect of the EW contribution is below 
1% of the tree-level QCD contribution. In the high tan f3 region the leading contributions 
come from the gg channel and are only partly canceled by the other channels. The EW 
contributions are of the order of a few per cents. In the ?7i S q Uar k = 900 GeV case partial 
cancellations between the 57 and the qq channels further suppress the EW contributions. 
In this case they are below 1%. The EW contributions to 6262 production are depicted in 
the right panels of Figure || and exhibit similar features. The NLO EW contributions are 
more pronounced since the corrections from the gg channel are more important. They are 
of the order of several percents, e.g. 5% for m squar k = 600 GeV and tan/3 > 45. 

It is worth noticing that the only practical effect of the resummation in the b/b sector 
is to change the value of the mass of the bottom squark. Indeed the impact of the effective 
Yukawa couplings in the computation of the amplitudes is negligible. This can be inferred 
from Figure 0, where we plot 



Aogw _ A(7 ew 



S a = eff n °- eff (A Q) 



'eff 

in the case of 62^2 production. Acr^ and Acr^^g denote the EW contributions with and 
without the effective Yukawa couplings, respectively. The overall effect of the effective 
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Figure 4: Relative difference between the EW contribution to the cross section for 6262 production 
with and without enhanced Yukawa couplings (black lines). The individual effect on the gg (66) 
channel is depicted by the rcd(-dottcd) lines. The qq channel is not affected by the resummation 
and hence not shown. ?7i squark is given in GeV. 
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Yukawa coupling, given by the black line, is in most regions below the permille level. The 
positive peaks for m, squar k = 600 GeV and m squar k = 900 GeV correspond to the zeros of 
the denominator in the definition ( [4,6[ ) , cf . Figure |3| 

4.4 Differential distributions 

Even though the EW contributions seem to have a rather small impact on the total cross 
section, they can become important in specific phase space regions. In Figure [| and |6| we 
consider differential distributions with respect to three different kinematical variables in 
the case of b\b\ and 62^2 production, respectively. We focus on the the SPSla' scenario 
and on the 14 TeV LHC. The left panels show to total EW contributions to the differential 
cross section. The tree-level EW contribution and the NLO EW contribution of the various 
production channels are depicted as well. The right panels show the impact of the EW 
contributions relative to the tree-level QCD cross section for each production channel. In 
contrast to the left panels, in the right panels the qq contribution is given by the sum of 
the tree- level and NLO EW contributions. 

Figure ||(a) refers to the transverse momentum distribution of the sbottom with highest 
Pt- Close to the threshold, i.e. in the region px < 300 GeV, the contribution of the 
gg channel is positive. Far from the threshold, this contribution becomes negative and 
relatively more important. In the low pt region the two contributions from the qq channel 
are different in sign, and their partial cancellation reduces the overall effect of this channel. 
In the high px region the qq channel increases its importance. The photon induced channel 
peaks at low pt and it is almost proportional to the LO QCD cross section, i.e. its relative 
yield is constant in px- As expected, the bb channel is irrelevant in the whole pt region. 
The total EW contributions have a small positive yield of the order of 1 — 2% in the low px 
region, while for px > 500 GeV the cross section is altered by 5 — 10%. It is interesting to 
note that a lower cut pxmin on the transverse momentum can significantly rise the impact 
of the EW contributions. For instance the cut pxmin = 320 GeV would discard the positive 
yield of the gg channel in the low px region. As a consequence the relative yield of the EW 
contribution to the total cross section would become of the order of —3.2%. 

The invariant mass distribution is displayed in Figure ||(b). The EW contributions 
exhibit the same high-energy behavior they have in the case of the px distribution. In 
this energy region they alter the leading-order prediction up to 10%. The peaks in the 
gg channel correspond to 62 and ti thresholds. Although in the low invariant mass re- 
gion no remarkable cancellations occur among the various channels, the overall positive 
contribution is small, of the order of 2%. 

In order to study how the EW contributions to the total cross section are altered by 
a lower cut Mi nV)TO i n on the invariant mass, we consider cr(Mi nVim i n ) defined as the total 
cross section integrated from the value Mj nv . m i n of the invariant mass. The upper-left 
part of Figure shows the relative yield of the total EW contributions to cr(Mi nVim i n ), 
together with the breakdown into the individual channels. The lower cut -Mi nVim i n excludes 
the region where the EW contributions are positive. Therefore the EW contributions 
decreases as Mi nV]m ; n increases, while their relative impact increases. For instance, in the 
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region M; nV]m i n > 1500 GeV the relative yield of the EW contributions <5 EW exceeds —5%. 
In this region, <5 EW is amplified by a factor of seven with respect to its value in the case of 
the fully inclusive cross section, cf. Table |3[ The cross section is reduced by a factor of five 
for Minv min = 1500 GeV (upper-right panel of Figure |?]). 

Figure |5](c) shows the pseudo-rapidity distribution, where always the squark with the 
higher absolute value of the pseudo-rapidity r\ is considered. The gap for zero rapidity is a 
consequence of this definition. The NLO EW contributions peak at \r)\ = 1 and dominate 
the EW contribution at this region. The contribution is negative for small values of \r]\ 
The total effect on the LO QCD cross section is up to 2%. 

Figure |6] shows the differential distributions for b^)\ production in the transverse mo- 
mentum (a), in the invariant mass (b), and in the pseudo-rapidity (c). In contrast to b\b\ 
production, the threshold behavior initiated by t<i is mild and hardly visible in both, the 
transverse momentum and the invariant mass distributions. The EW contribution is small 
and its relative yield stays below 5%, even in the high energy region. This is expected since 
in the SPSla' scenario 62 is mostly right-handed. Interestingly, the contributions from the 
qq channel and the 57 channel almost cancel in most parts of the phase-space. Therefore 
the EW contributions are well approximated by the gg channel corrections. The lower 
panels of Figure @ show 0"(Mm Vjm i n ) in the case of 62^2 production. In this scenario the EW 
contributions to 62^2 production can be safely neglected for each value of Mi nVjm i n . Even 
in the case Mi nVim ; n = 2 TeV the EW contributions change the Born cross section only by 
an amount of the order of 2%. 

In Figure ^ we consider the differential transverse momentum distributions for b\b\ 
production in the SPS8 scenario for 14 TeV (upper panels) and 7 TeV (lower panels). It 
is worth analyzing the EW corrections to the gg channel in this scenario. The sbottom 
mass is heavier in the SPS8 scenario than in the SPSla' scenario. Therefore the gluons 
producing a sbottom-anti-sbottom pair have a bigger typical momentum fraction x in the 
SPS8 scenario than in the SPSla' scenario. Since the gluon PDF falls off rapidly at high x 
values, the negative EW contributions of the gg channel in the high px region are strongly 
suppressed in the SPS8 scenario. In the VS = 7 TeV case the typical value of x it is even 
bigger and this phenomenon is enhanced. In the SPS8 scenario the cancellation between 
the positive low px corrections and the negative high pr ones is less effective. Therefore, 



as mentioned in Section 4^, the overall (positive) gg channel contributions get relatively 



enhanced in the ^/S = 7 TeV case. 
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Figure 5: Differential distributions for the transverse momentum pr, the invariant mass M; nv and 
the pseudo-rapidity r\ for b\b\ production at the 14 TeV LHC within the SPSla' scenario. Shown 
are the tree-level and NLO EW cross section contributions for the various production channels (left) 
and the impact of the NLO EW contributions relative to the LO QCD cross section (right). In the 
left panels the tree-level and NLO EW contributions for the qq channel are plotted separately. In 
the right panels they are treated inclusively. 
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Figure 6: Same as Figure |g but considering 62 &2 production. 
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Figure 7: Left: Relative yield of the EW contributions of the different production channels in 
cr(Minv,min)- Right: LO QCD prediction for the same observable. 
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Figure 8: Differential transverse momentum distribution for bib\ production within the SPS8 
scenario for the 14 TeV (upper plots) and 7 TeV (lower plots) LHC. 
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Figure 9: Relative yield of the various liadronic processes with a bottom-squark pair in the final 
state as a function of the common squark mass breaking parameter m squark . The parameter m squar k 



is defined in Section 4.3 



5. Other processes leading to bottom-squark pair production 

In this section we describe the remaining processes of (JO]) leading to a pair of bottom- 
squarks in the final state. The processes are non-diagonal sbottom-anti-sbottom produc- 
tion and (anti-)sbottom-(anti-)sbottom production, 

PP^kblhbl; PP^b a b^b* a %- a,/3e{l,2}. (5.1) 

The leading-order cross sections, of the order 0(a^), are given by 

r<n+ f 
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where ro(a,/3) = (m^ + ) 2 /<S is the production threshold. The parton luminosities 



dLy/dr are defined according to Eq. ( |2.1| ). The tree- level EW contributions are of the 
order 0(a s a) and 0(a 2 ) and read as follows, 
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The numerical impact of these processes on bottom-squark pair production at tree-level 
is rather small. In Figure ^ we show the relative yield of diagonal sbottom-anti-sbottom 
production, Eq. (|l.2| ), and of all the processes described in this section as a function of the 
soft breaking parameter m squark . Owing to the smallness of the (anti-)bottom PDF, the 
tree- level contribution of the processes (|5.1| ) is below 1%. We therefore do not include NLO 
EW corrections; they are expected to be even smaller and thus can be safely neglected. 
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6. Conclusions 

We have studied the EW contributions to 6-squark pair production at the LHC within the 
MSSM. The tree-level EW contributions from both the gg-annihilation and the gluon fusion 
channels have been supplemented by the contribution from the photon-induced channel, 
not included in previous discussions. We have presented the first complete computation of 
the NLO EW contributions to diagonal sbottom-anti-sbottom production; together with 
the QCD corrections they complete the NLO analysis of the class of squark-antisquark 
production processes. 

Renormalization of the 6-sector is shown in detail, with a check of the reliability of the 
adopted renormalization scheme in the numerical analysis. The potentially large corrections 
to the bottom Yukawa couplings have been resummed, introducing effective couplings. In 
the scenarios considered, the main effect is to change the value of the masses of the bottom 
squarks. 

The EW contributions to the total cross section are strongly scenario dependent. How- 
ever, in all the scenarios considered they are of the order of few percents of the LO con- 
tribution in inclusive cross sections. Their size is partly due to strong cancellations among 
different channels. The EW contributions of different channels peak in different regions of 
the phase space. Therefore, the EW contributions are enhanced whenever kinematical cuts 
are applied. For similar reasons the impact of the EW contributions is more important in 
differential distributions, in particular in the high-energy region. In the SPSla' scenario, 
and in the case of b\b\ production, they can even exceed 10% of the LO contributions for 
invariant mass and transverse momentum distributions. 
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A. Renormalization of the stop and sbottom sector 



In the MSSM the kinetic terms in the stop and sbottom sector read as follows 

o^qR ) \qR 



£=E Wl,9^r) ( ) ~ (£,ft)Mf ( f ) I (A.l) 



q=t,b v 

with the squared-mass matrix 

M 2 = (Ml + m\ + m| cos 2/3 (T ? 3 - e q s 2 w ) m q (Aq - fiXg) \ 

5 ^ m g (^ ? --yuA 9 -) M| fl + m2 + e 9 Mf cos 2^4 J ' l " J 

Mi, Mq t ft, Ag are soft breaking parameters, while is the super symmetric Higgs mass 
parameter. m q , e q , T 3 are the mass, the charge and the isospin of the quark q, respectively. 
s w is the sine of the weak mixing angle, while A^ = cot /3 and A^ = tan /3. The matrix (A. 2) 
is symmetric and can be diagonalised by an orthogonal matrix H q such that 

(«A=rJ« l ), K q = ( R ^ iR ^A. 

In the rotated basis the squared-mass matrix is diagonal 

The rotation matrix ~R q can be parametrized in terms of a mixing angle 6 q 

/co S «, si„«-\ or 
y — sin o q cos v q J 

depending on the sign of the determinant of R^. The mixing angle 6 q and the trilinear 
coupling A q are related via 

sin & = 2m ^ ~ ^ e e = det [R,] . (A.3) 

Because of SU(2)-invariance the eigenvalues of M| and M~ are connected. In particular 
they have to satisfy the relation 

^i,i m L + R ~b2,i m l 2 - m b = R \\,\ m \ + R i2,i m h - m t- m w cos 2/3. (A.4) 



Definition of the "DR bottom-quark mass" renormalization scheme 

In this scheme the independent parameters are chosen to be 

,2 2 2 



i> m t 2 ' "V m *' mb ' ^6' ^ ( A ' 5 ) 
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The squark mass squared and the top quark mass are defined in the on-shell (OS) scheme, 
while the bottom quark mass and the trilinear coupling A b are fixed using the DR pre- 
scription, 
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<5/i and <5tan/3 appearing in Eq. ( A.6 ) are defined in the DR scheme. S dlv ' is the divergent 
part of the scalar self energies defined according to the following Lorentz decomposition, 



The stop mixing angle is defined according to 
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In the "DR bottom-quark mass" scheme Ai, 6r and are dependent quantities. Their 



counterterms read as follows 
1 



8Ar = — 



R tl,l R tl,2 ( 5m ~ tl ~ 5m t 2 ) + 6 ( R tl,l R t2,2 ~ R tl,2 R t2,l) 50 t( m t 1 



rrtr 



-{Ai — //cot (3)5m t 



+ 5 /j, cot f3 — //cot /3<5tan/3, 



(A.9) 



61,1 



—2m t 5m t — cos 2/3 — m^/<5 cos 2/3 J + i?^ 2 5m?^ 
+2R~ bl i R ~ b i 2 ~~ 5// tan /3 — //J tan /3) 
+5m b (2R~ bhl R ila (A~ b - /ztan/3) + 2(1 - 2R\ ia )m h 



(A.10) 



R bl,l R b2,2 + R bl,2 R b2,l 

+m b {5A~ b — fi5 tan /3 — 5[i tan /3) 



R bi,i R bi,2 { dm l 2 ~ 6m bl ) + 0*5 ~ [i tan P)6m b 



(A.ll) 
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The mass of the W boson, m\y, is renormalized on-shell. In order to get finite Green 
functions, we also need the wavefunction renormalization of the bottom and sbottom, 



Q2 



bare 
bare 



u±q 



bare 




6Z S 



(A.12) 



(A.13) 



In the case of bb* production only the diagonal entries of the matrix 5Zq are needed. They 
are defined as 

d£(fc 2 ) 
dk 2 Ifc2=m? 



8Z, 



qa,a 



-Re{£^ Q (m 2 )}, S'(m 2 ) 



The wavefunction renormalization constants of the bottom-quark read as follows 

5Z bL/R = -Re {z bL/R (m 2 ) + m 2 (v' bL {m 2 ) + Z' bR (m 2 b ) + 2£' bS (m 2 )) } . (A.14) 

In the processes considered, b\ is an external particle and its mass has to be defined 
on-shell. Therefore we set the value of the 61 mass to its OS value, obtained using the 
following relation 

m L os = m l + H " Re { E &i,iK)} • (A. 15) 

The renormalization constant of m~ enters the counter terms in the last diagram of Fig- 



ure B.3 and it is fixed in accordance with our choice of the 61 mass, 



5m\ 



1, OS 



Re { s Si,iH,os)} 



(A.16) 



The renormalization constants have to be evaluated at 0{a) and enter the calcula- 



tion via the counterters in Figure |B.3| and Figure |B.4j . The explicit expressions for the 
counterterms are given by 



b a 



P 



ba,a 



5m 2 




b a 



ig 2 (T^+T^)5Z- baagflu , 



9 

JLfiJUULSJix " 



-i9sT c 5Zr(k + k')^, 




-ig s T c (5Z qL -f^uj- + SZqRip 
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B. Feynman diagrams 



In this Appendix we list the Feynman diagrams relevant for b a b* a production at tree-level 
and at next-to leading order electroweak. S (S^) denotes the neutral (charged) scalar 
Higgs and Goldstone bosons. V = 7, Z. 




Figure B.l: Tree-level Feynman diagrams for b a b* a production, (a) and (b) show the QCD dia- 
grams for the gg, qq and bb channels, (c) are the EW diagrams. They are not present in the gg 
channel at tree-level, (d) is the E W tree-level diagram involving CKM matrix. It is only present 
for initial uu or cc quarks. 




K K * b a \ X l a 



7 



7 



Figure B.2: Lowest-order Feynman diagrams for the gluon-photon fusion process 57 — > b a b* a 
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Figure B.3: Virtual corrections to the process gg —> b a b* a . A common label V is used for the 
neutral gauge bosons "f,Z°, while S denotes any of the neutral Higgs bosons or the neural Goldstone 
boson h° 7 H° 7 G , and denotes the charged ones i? ± ,G ± . Crossed diagrams are not shown 
explicitly. The diagrams containing the countcrtcrms arc depicted in the last line. The counterterms 
have to be evaluated at 0(a). 
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q 9 





Figure B.4: Virtual contributions to the process qq — > b a b a . The diagrams result from EW 
insertions to tree-level QCD diagrams and from QCD insertions to tree-level EW diagrams. V, S, 



and are defined as in Figure B.3. Crossed diagrams are not shown. The counterterms in the 
last two diagrams have to be evaluated at 0(a). 





Figure B.6: Real photon emission. The first six diagrams correspond to the gg channel while the 
last three correspond to the qq channel. 
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(a) 




(I)) 

Figure B.7: Real gluon emission for the qq channel, (a) QCD based diagrams, (b) EW based 
diagrams. 




Figure B.8: Feynman diagrams contributing to real quark radiation, (a) QCD based diagrams, 
(b) EW based diagrams. 
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